Microscopic formula for transport coefficients of causal hydrodynamics 



T. Koide 1 

1 Instituto de Fisica, Universidade Federal do Rio de Janeiro, C. P. 68528, 21945-970, Rio de Janeiro, Brazil 

The Green-Kubo-Nakano formula should be modified in relativistic hydrodynamics, because of 
the problem of acausality and the breaking of sum rules. In this Rapid Communication, we propose 
a formula to calculate the transport coefficients of causal hydrodynamics based on the projection 
operator method. As concrete examples, we derive the expressions for the diffusion coefficient, the 
shear viscosity coefficient, and corresponding relaxation times. 
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The hydro-dynamic model has been applied to analyze 
collective aspects of relativistic heavy-ion collisions. The 
analysis, so far, has been implemented extensively for 
ideal fluid. However, there is no theoretical reason that 
we can ignore the effect of dissipation. The relativistic ex- 
tension of dissipative hydrodynamics is not trivial [3, 0] • 
Because of the problem of acausality and the breaking of 
sum rules, the simple covariant extension of the Navier- 
Stokes equation does not work. To obtain a consistent 
theory, we need to introduce a memory effect with a re- 
laxation time [1,01, HI- Then, the diffusion-type equation 
is modified to the telegraph-type equation. 

This modification also affects the definition of the 
transport coefficients. In the usual non-relativistic hydro- 
dynamics, it is known that such transport coefficients can 
be calculated by using the Green-Kubo-Nakano (GKN) 
formula However, in relativistic hydrodynamics, 

there is a problem of using the GKN formula to esti- 
mate the transport coefficients of causal hydrodynamics. 
For example, let us consider the diffusion process which 
is described by 



d 

—n(x,t) = D ac V 2 n(x,t). 
at 



(1) 



Then, by using the GKN formula, the diffusion constant 
is given by || 

Dac = (?£) lim lim [°° dt [ d 3 xe~ 1 ^-^ 
\dnj ^-om J J 

1 f - 

x-J d\{j(x,t-i\)j(O,0)) eq , (2) 

where \i represents a chemical potential and ( ) eq denotes 
the equilibrium expectation value with the density matrix 
p = e- f)H /Z using Z = Tr[e- 0H ]. Here, j(x,t) is a 
current operator associated with n(x, f) (see Eq. 
below) . In deriving this expression, we assumed that the 
coarse-grained dynamics of n(x, t) is given by Eq. (P). If 
the coarse-grained dynamics is changed, we cannot use 
Eq. ©. Thus the GKN formula should be modified. 

As a matter of fact, it is known that the diffusion equa- 
tion does not obey relativistic causality: the propaga- 
tion speed of signal of Eq. ([I]) exceeds the speed of light 
[1, 01 IH- Moreover, the diffusion equation is inconsistent 
with sum rules associated with conservation laws [f| 0] • 



Aiming to avoid these problems, a frequently considered 
alternative is given by the telegraph equation, 

J^n(x,i) + ~n(x,t) = D c \7 2 n(x,t). (3) 

Here, the parameter r is the relaxation time and it is 
necessary to obtain a causal theory which is consistent 
with sum rules [y, |7| . When the coarse-grained dynamics 
is given by Eq. ([3]) instead of Eq. (Q}, the causal diffusion 
coefficient D c is not given by Eq. @ . It should be noted 
that the telegraph equation will not be a unique solution 
to solve the problem of acausality. See § for details. 

In Rapid Communication, we derive a general formula 
of transport coefficients of causal hydrodynamics. For 
example, as we will see later, in the case of the causal 
diffusion process described by Eq. ([3]), the causal diffu- 
sion coefficient and the relaxation time are given by Eqs. 
(fTT)) and (HHJ). Interestingly, the diffusion constant ([2]) 
does not represent the causal diffusion constant. It re- 
duces to the inverse of the relaxation time under some 
approximation. 

First of all, we consider an operator of a gross (hydro- 
dynamic) variable O(x). The time-evolution of the O is 
given by the Heisenberg equation of motion, 



d_ 

df 



0(x,t) = iLO(x.,t), 



(4) 



where L is the Liouville operator defined by LA = [H, A] . 
Here, A is an arbitrary operator and H is the Hamilto- 
nian of our system. The operator O is arbitrary, but 
should satisfy {0} eq = @, @. The Heisenberg equa- 
tion of motion contains the information not only of gross 
variables associated with hydrodynamic time scales, but 
also of microscopic variables with shorter time scales. To 
carry out coarse grainings of the latter irrelevant infor- 
mation, we introduce a projection operator P and its 
complementary Q = l-P 0, 1 E, El B O • Then, 
the Heisenberg equation of motion is rewritten as 

^0(x,t) = e lLt PiLO{i) 

+ f dse tL{t - s) PiLQe lLQ HLO{-x)+Qe lLQt iLO{-K).{h) 
Jo 

This equation is called the time-convolution equation 
[H, EH- The first term on the right-hand side (rhs) is 
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called the streaming term and corresponds to a collective 
oscillation which does not produce entropy. The second 
term represents the dissipative part of the time evolu- 
tion. The third term is interpreted as the noise term. 
The irreversible current is described by the second term. 
Irreversible currents are proportional to thermodynamic 
forces and the transport coefficients are defined by the 
proportional coefficients. Thus the first and third terms, 
which are irrelevant to the definition of transport coeffi- 
cients, are ignored in the following discussion. 

So far, we did not specify the form of the projection op- 
erator. In this work, we use the Mori projection operator 



PA 



d 3 xd 3 x'(A, O t (x))(0(x), 0t (x'))~ 1 0(x'), (6) 



where the Kubo's canonical correlation is defined by 



Jo P 



eq ■ 



(J) 



For the physical meaning of the Mori projection, see 0, 
ITTj |. Substituting the Mori projection into Eq. ([5]), we 
have 

ft 



d_ 

df 



0(x,i) = -^ ds J dVS(x,x';s)0(x',t-s). (8) 
The memory function 2 is given by 
S(x,x";t) = -6{t) J d 3 x' 

x(iLQe 4iQt iLO(x),O t (x'))(0(x'),O t (x"))" 1 . (9) 

It is difficult to calculate the memory function because 
of the coarse-grained time-evolution operator e ji(3t . Usu- 
ally, we approximately replace this term with the usual 
time-evolution operator e lLt . As is discussed in [l5j], the 
approximation was considered to be justified at least for 
the low momentum limit, and the diffusion- type equation 
is derived. However, as is shown in [7|, this approxima- 
tion is not justified even for the low momentum limit, and 
gives rise to several problems: inconsistency of causality 
and sum rules. We will come back to this point later. 

Recently, the coarse-grained time-evolution operator 
was calculated exactly @, H, E3|- By using the same 
method, we find that the Fourier transform of the mem- 
ory function is given by 



Z F (k,u) = - 



X F (k,co) 
l + X F (k,u) 



X F (k,t = 0)X F (k,u 

l+X F {k,Uj) 



{.10) 



Here, X F (k, ui) and X F (k, w) are given by the Fourier 
transforms of the following correlation functions, respec- 
tively, 



X(x-x",t) 
X(x-x",t) 



--0(t) J dV<[0(x,t),Ot (x '])) e(J 
x(0(x'),O t (x"))- 1 , 



-~9(t) J d 3 x'([iLO{x,t),0\x!)]) eq 

x^x^ovor 1 - 



(ii) 

)eq 

(12) 



This is the first main result of this work. A similar expres- 
sion of the memory function is obtained for the current- 
current correlation function in [3], although the term 
that corresponds to the second term on the rhs of Eq. 
(|T0| does not exist because they considered the case of 
X F (k,t = 0) = 0. 

To obtain the expression of transport coefficients, 
we have to employ an approximation to derive a 
time-convolutionless equation from the integrodifferen- 
tial equation. Following @, Q , we separate the memory 
function into two terms. 



0_ 

Of 



0(k,t) = - [ drf7 2 (k,<-r)0(k,r) 
Jo 

- I dr<f>{k,t-T)0{k,T), (13) 
Jo 



where the frequency function and the renormalized mem- 
ory function are, respectively, defined by 



n 2 (k,t) = l J ^Im[H F (k,o;)] e 
*(M) = J ^Re[S^(k, w )]e- 



-iurt ( M ) 
luJt . (15) 



To employ this approximation, we have to know the 
temporal behavior of the two functions. So far, the be- 
haviors of the two functions are investigated for the chiral 
orderparameter Q and the nonrelativistic number den- 
sity [3]. For these cases, the two functions showed com- 
mon properties; the frequency function converges to a 
finite value and the renormalized memory function van- 
ishes at late time. Thus we employ an important as- 
sumption for the temporal behavior of the two functions. 
The renormalized memory function relaxes rapidly and 
vanishes at late time. On the other hand, the frequency 
function converges to a finite value after short time evo- 
lution. Then, the approximation means to replace the 
frequency function and the renormalized memory func- 
tion with the following functions: 



n 2 (k,t) -^D k k 2 , 



Tk 



where 



1 



-5 lim fi 2 (k,i), 



dt$(k,i). 



(16) 

(17) 
(18) 



These two definitions are the second main result of this 
Rapid Communication. The coefficient £>k gives us the 
dissipative coefficient and the coefficient Tk is the relax- 
ation time. As a matter of fact, by using these expres- 
sions, we arrive at 



d_ 

at 



i /•* 

Oik, t) + —0(k, t) = -L> k k 2 / drO{k, r). (19) 

Tk Jo 

This equation has the same form as the causal diffusion 
equation ([3]) except for the momentum dependence of the 
relaxation time Tk Q ■ 



3 



The same equation as Eq. (QjJJ has already been given 
by 0)01 • However, the perturbative approximation was 
applied in the calculation of the memory functions there. 
The result of this Rapid Communication is exact except 
for the local approximation. 

The validity of the above assumption for the frequency 
and renormalized memory functions should be checked by 
calculating these temporal behaviors, numerically. How- 
ever, we can see that the assumption is reasonable by 
usingthe final value theorem of the Laplace transforma- 
tion [Tj, [9j . When the renormalized memory function con- 
verges to zero at late time, its Laplace transform should 
satisfy, 



lim $(x, t) = lim s$ /j (x, s) = 0. 

t— >oo a— >0 



(20) 



Similarly, the dissipative coefficient which is calculated 
from the frequency function should be given by 



Duk 2 



lim n 2 (k,t) 



lims(f2 2 ) L (k,s). (21) 

s— »0 



Here, the suffix L means the Laplace transform. 

To obtain the dissipative coefficient and the relaxation 
time, we have to know the correlation functions X F (k, ui) 
and X F (k, us). Both of them can be calculated by us- 
ing the imaginary-time Green function, C/(x, t;x',t') = 
— Tr{pT T (3(x, r)0' l '(x, r')}. Once we have it, the correla- 
tion functions X F (k, u>), which is essentially the retarded 
Green function, can be estimated by using the analytic 
continuation of the imaginary-time Green function. The 
other correlation function X F (k, u) can be calculated by 
using the following relation,: 

J^X(x - x"; t) = l(x - x"; f) + Z(x - x"; t)S(t). (22) 

To see the relation between our formula and the GKN 
formula, we would like to discuss the diffusion process 
of a conserved number density n(x, t), by setting O = 
<5n(x) = ri(x) — (n(x)) eq . When we approximately replace 
e tLQt Q f q witfo e iLt ag wag men tioned below Eq. 

([5]), the memory function is given by 



(k,w) w -X*(k,u). 



(23) 



Here, we used X F (k,t = 0) = for O = <5n(x). The 
correlation function X F (k,t) is 

x F (k,t) = <w(^) 

x J d 3 xe- ikx ifc 2 (J(x,i),J(O,0)), (24) 

where p is the chemical potential. Here, we intro- 
duced the current operator by dn(x.,t)/dt = — Vj(x, t). 
And the equal-time correlation function is approximately 

given by / d 3 xe-^(5n(x, 0), Sn(0, 0)) « \ (fa) @. 
At the low momentum limit, we have 



D k = 0, — - D ac k 2 , 

Tk 



(25) 



where D ac is defined by Eq. ([2]). Then, one can see that 
Eq. (Til)|) is reduced to the diffusion equation ([T]) , and the 
diffusion constant is given by Eq. ((2]). From this anal- 
ysis, we found that our formula is reduced to the GKN 
formula when we ignore the effect of coarse grainings on 
correlation functions. Then, the GKN formula gives not 
the diffusion coefficient D c but the inverse of the relax- 
ation time l/rfc. 

It should be noted that the diffusion equation cannot 
satisfy the f-sum rule. In general, when we take the ef- 
fect of the coarse-grained time-evolution operator into 
account, £>k does not vanish and we obtain the telegraph 
equation. The telegraph equation can be consistent with 
the f-sum rule 0, 0| . 

Next, we apply our formula to a nonrelativistic fluid 
to derive the shear viscosity constant. In the classical 
calculation, the commutator should be replaced by the 
Poisson bracket. For simplicity, we consider a particu- 
lar case of shear flow, where the fluid velocity u x (y,t) 
points in the x direction and varies spatially only in 
the y direction. Then, the current density is given by 
J x {hy) = P u x = J2pfe' lkyVi , where p is a mass density 
and pf is a momentum of a constituent particle. By set- 
ting O = J x (k y ) and substituting into Eq. (fl9)) . we have 
the evolution equation of the current density, 



lr(k v ,t) + —r(k y ,t) = 

0t T k 



-D k „ky I dTpu x {k Vl - 
o 



(26) 

If we ignore the effect of the coarse graining, we obtain, 

1 B f°° 

D k = 0, - = k 2 ^- / dt(P xy {k,t)P xy {~k, 0)) eq ,(27) 
Tk pv Jo 

where the pressure tensor is defined by J x = ik y P xy (k y ) 
and V is the total volume of the system. At the low 
momentum limit, the inverse of the relaxation time is 
reduced to the GKN formula for the shear viscosity con- 
stant [11]. 

Similarly, we can derive the formula for the shear vis- 
cosity of the relativistic dissipative hydrodynamics 0. 
We consider a particular situation which is the same 
as the nonrelativistic case discussed above. Then, the 
energy-momentum tensor obeys the following equa- 
tion of motion: ^T 0a; = —§^Ty x = -^ xy , where 
is the shear viscosity. Here, we ignore the reversible part. 
As is discussed in 0] , the shear viscosity is given by solv- 
ing the telegraph-type equation. Then, the linearized 
equation of motion in terms of the fluid velocity is 



o 2 



T 



0x 



1 d 



T 



Or 



d 2 



o. 



dt 2 rdt 2rdy 2 

On the other hand, by setting O = T 0x , we obtain 

^-T 0x (k y ,t) + ±T 0x (k y ,t) 



(28) 



= —D 



k 2 



dTwu x {k y ,T), (29) 
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where w is an enthalpy density. Here, we assume that the 
fluid velocity is given by u x (k y ) — T 0x (k y )/w in the lin- 
ear approximation of the fluid velocity. Comparing with 
Eq. (f2"5| . we obtain the expressions of the transport co- 
efficients of relativistic dissipative hydrodynamics. The 
relaxation time and the shear viscosity constant are given 
by 

7] = lim 2wD k r k , (31) 

respectively. In particular, the correlation function 
X(k, t) is given by 

X F (k,t) = 9(t)^ y*d 3 a;d 3 a;ie- lkx fc 2 (T Q/3 (x,i),T Q/ 3(x 1 ,0)) 

xCT -^, 0)^(0, 0))- 1 (a, (3= 1,2,3). (32) 

If we ignore the effect of the coarse-grainings again, our 
calculation reproduces the shear viscosity given by [i~7l ]. 
In this sense, our expression of the shear viscosity is the 
generalization of the well-known result. 

In this work, we derived a formula to calculate trans- 
port coefficients of causal hydrodynamics. In our for- 
mula, we took the effect of the coarse grainings in the 
memory function into account. Then, the coarse-grained 
equation is given not by the diffusion-type equation, but 
by the telegraph-type equation. By comparing our equa- 
tion with phenomenological equations, we obtained the 



expressions of the diffusion coefficient, the shear viscosity 
coefficient and the relaxation times. This formula repro- 
duces the GKN formula when the effect of the coarse 
grainings is ignored in the memory function. 

Near phase transitions, the relaxation of the fluctua- 
tions of thermodynamic quantities is decelerated. This is 
called the critical slowing down and is expected to hap- 
pen for the shear viscosity [HI . We have already showed 
that the critical slowing down near the chiral phase tran- 
sition can be described in this formalism 0. 

In this Rapid Communication, we derived the linear 
equation 1(5)1 and ignored the effect of nonlinearity in- 
cluded in the noise term, the third term of Eq. ([5]) , which 
is negligibly small near equilibrium. When we discuss the 
far-from-equilibrium dynamics, which will be described 
by the nonlinear hydrodynamics, we have to extend the 
definition of the projection operator so as to extract the 
non-linear effect from the noise term 0, [ll[ . 

The formula of transport coefficients of extended ther- 
modynamics is discussed also by Ichiyanagi [19]. The 
formula, however, contains the coarse-grained time- 
evolution operator, explicitly. Thus to apply the formula 
to practical problems, we need to rewrite the expression 
as was done in this work. 

The author thanks G. S. Denicol, F. Gelis, T. Kodama, 
and Ph. Mota for helpful discussions. This work was 
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